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Abstract 

We consider the following problem: find a fixed-order linear controller that maxi- 
mizes the closed-loop asymptotic decay rate for the classical two-mass-spring system. 
This can be formulated as the problem of minimizing the abscissa (maximum of the 
real parts of the roots) of a polynomial whose coefficients depend linearly on the 
controller parameters. We show that the only order for which there is a non-trivial 
solution is 2. In this case, we derive a controller that we prove locally maximizes 
the asymptotic decay rate, using recently developed techniques from nonsmooth 
analysis. 



1 Problem Statement 



We consider the system shown in Figure Q consisting of two masses interconnected by a 
spring, a typical control benchmark problem which is a generic model of a system with 
a rigid body mode and one vibration mode [10 . If the first mass is pulled sufficiently 
far apart from the second mass and suddenly dropped, then the two masses will oscillate 
until they reach their equilibrium position. 



The control problem we study in this note consists of appropriately moving the second 
mass so that the first mass settles down to its final position as fast as possible; more 
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Figure 1: Two-mass-spring system. 



specifically, we want to maximize the asymptotic decay rate. For this we use a linear 
feedback controller between the system output (measured position of the second mass) 
and the system input (actuator positioning the first mass). This control problem can be 
formulated as the minimization of the abscissa (maximum of the real parts of the roots) of 
a polynomial in the complex Laplace indeterminate, whose coefficients depend affinely on 
the controller parameters. This polynomial is the denominator of the closed-loop system 
transfer function. 

For notational simplicity, we assume here that both mass weights and the spring constant 
are normalized to one. As shown in ^H], the polynomial is then given by 

P(a?, V)( s ) = ( s4 + 2s 2 )a;(s) + y(s) 

where x(s) and y(s) are respectively the denominator and numerator polynomial of the 
controller transfer function to be determined. This transfer function is assumed to be 
proper, i.e., 

m = deg x(s) > deg y(s). 

The integer m is called the order of the controller. Without loss of generality, we take x(s) 
to be monic. Letting P" denote the linear space of polynomials with complex 1 coefficients 
and degree < n, we therefore write y(s) G P' m and x(s) — s m G P m_1 . 

For a given polynomial q G P n , we define the abscissa of q by 

a(q) = max {Re z : q(z) = 0} . 

Since the roots of p(x, y) are the closed-loop system poles, the closed-loop two-mass- 
spring system is stable if and only if a(p(x, yj) < 0. We are interested in maximizing the 
asymptotic decay rate of the system, i.e., solving the optimization problem 

inf a(p(x,y)). (1) 

y(s)SP m , i(s)-s m 6P m " 1 

Section |21 shows that for order m < 1, there is no stabilizing controller, i.e., no (x,y) 
such that a(p(x,y)) < 0. Section El shows that when m > 3, the abscissa a(p(x,y)) is 
unbounded below. Section 0] studies the more interesting case m = 2, and gives a formula 
for (x, y) that, in Sectional we prove is a strong local minimizer of a(p(x,y)). Sectional 
plots the time response of the optimized controller and discusses the issue of robustness. 
Concluding remarks are made in Section [7| 

work with the space of polynomials with complex coefficients for technical reasons; the chain rule 
that we use in Section [S] is most naturally stated in this context. 
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2 First-Order Controller Design 



If we assume that the controller has order m — 1, then both x(s) and y(s) are first degree 
polynomials, say 

x(s) = x + s, y(s) =y + yis, 

with 

p(x, y)(s) = y + yts + 2x s 2 + 2s 3 + x s 4 + s 5 . 

This polynomial is stable, i.e., it has all its roots in the open left half-plane, if and only 
if all the principal minors of its Hurwitz matrix 

x Q 1 

2x 2 x 1 

y y 1 2x 2 x 

y yi 2x 

. y _ 

are all strictly positive, see e.g. [Hj. This can never be the case since the 2-by-2 northwest 
minor has rank one for all Xq. Hence a controller of first order (or less) cannot stabilize 
the two-mass-spring system. 

Note that in problem the minimum abscissa a(p(x,y)) = is attained for any static 
controller x(s) = 1, y(s) = yo = k with k e [0, 1], because then 

p(x, y)(s) = k + 2s 2 + s 4 = (s 2 + 1 + y/l - k)(s 2 + 1 - y/l + k) 

has only imaginary roots. Similarly, a(p(x, yj) = is attained for any first-order controller 
x(s) = x Q + s, y(s) = yo + 2/i s such that xq = yo = 0, y\ = k and k e [0, 1], since then 

p(x,y)(s) = (k + 2s 2 + s 4 )s. 

3 Third-Order Controller Design 

If we seek a controller of order m = 3 then we can write 

6 

p(x,y)(s) = (s 4 + 2s 2 )(x + x x s + x 2 s 2 + s 3 ) + y + y x s + y 2 s 2 + y 3 s 3 = s 7 + ^2 / p i s l . 

i=0 

By identifying powers of the indeterminate s, we derive the linear system of equation 












1 








" 




" xo ' 




Po 














1 










Xi 




Pi 


2 














1 







x 2 




P2 





2 














1 




yo 




P3 


1 





2 
















2/i 




Pi 





1 



















2/2 












1 
















2/3 




PS 



3 



This 7-by-7 matrix is called the Sylvester matrix of polynomials s 4 + 2s 2 and 1, and it is 
non-singular because these two polynomials share no common roots. In other words, we 
can find controller coefficients defining any desired closed-loop polynomial. This is what 
Dorato jH] calls the fundamental theorem of feedback control, namely the fact that the 
poles of a single-input-single-output linear system of order n can be placed arbitrarily by 
a linear controller of order n — 1. Consequently, an arbitrarily large negative abscissa can 
be achieved in closed-loop by a controller of order three (or more). 

For example, by solving the linear system of equations shown above, we obtain that the 
controller polynomials 

x(s) = (-35z 3 + Uz) + (2lz 2 - 2)s - 7zs 2 + s 3 , 

y(s) = -z 7 + 7z 6 s + (-21z 5 + 70z 3 -28z)s 2 + (35z 4 - A2z 2 + 4)s 3 

place all the closed-loop poles at an arbitrary real value z. It follows that problem (JIJ is 
not bounded below for m > 3. 



4 Second-Order Controller Design 

In the case of a controller of order m = 2, we have 

p(x, y)(s) = 4 + 2s 2 )(x + xxs + s 2 ) + y + y x s + y 2 s 2 . (2) 

We can cluster all the closed-loop poles at a real negative value z by solving the following 
system of equations 

y = z 6 
yi = —6z 5 
y 2 + 2xq = 152 4 
2xi = ~20z 3 
2 + x = Ihz 2 
X\ = —6z. 

We observe that the only constraint on z is enforced by the fourth and sixth equations, 
namely 

5z 3 = 3z. 

We rule out the case z = since the closed-loop system would be only marginally stable, 
and we extract the negative solution 



1 ' -0.7746. 
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The controller coefficients can now be derived by substitution, resulting in 
* _ 7 * _ ^ » _ 27 „ _ 54y/l5 . _ ^43 

x -(, Xl - 5 , 2/o- 125 , Vl - 125> y 2 - 5 > W 
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yielding 



/ * *w n 6 6V15 5 n 4 12V15 3 27 2 54V15 27 
p(x*, y*)(s) = s 6 + — s 5 + 9s 4 + -^s 3 + -s 2 + -^-s + — . (4) 

We should note that we realized that all roots could be clustered at a single point z* after 
performing numerical experiments using hifoo jSl Section 6], a new toolbox for low order 
controller design using methods of nonsmooth optimization. 



5 Local Optimality Certificate 

In this section we prove that z* = — is, at least locally, the minimal abscissa that 
can be achieved with a second-order controller, i.e., a local minimizer of problem ([TJ for 
m — 2. This is nontrivial, since one might think it is necessary to consider all possible 
splittings of the multiple root under perturbation. We prove local optimality using recent 
advances in nonsmooth analysis. 

Recall that P n is the linear space of complex polynomials with degree less than or equal 
to n (with complex dimension n + 1) and let C n denote the space of complex vectors of 
length n. We write elements of C n as row vectors. 

We show that (x*,y*) defined by Q locally optimizes the abscissa of p(x, y), in the sense 
that any sufficiently small perturbation to Q strictly increases the maximum of the real 
parts of the roots. In fact, we prove that © is a sharp local minimizer, in the following 
sense. 

Theorem The abscissa of the polynomial p(x, y) defined in J||) is locally minimized by 
(x*, y*) with coefficients given by Furthermore, for (x, y) sufficiently close to (x*, y*), 
we have 

a(p(x, y)) > a(p(x*,y*)) + r\\d\\ 
where r is a positive constant and 

d=[x - x* ,x 1 - x\,y Q - 2/0,2/1 - V*i,V2 - y* 2 ] ■ 



The proof of this theorem is the subject of the rest of this section. It follows j3J Section III 
A] quite closely. We start by making a change of variables to the polynomial p(x,y)(s), 
namely 

t = s — z , 

Q = [<?o, qi] = [xq, Xi\ - [x* , x{], 

r = [r„, n, r 2 ] = [y , yi, y 2 ] - [y* , v\ vl\- 

A few lines of Maple [7j, specifically 
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p:=(s~4+2*s~2)*(s"2+xl*s+x0)+(y2*s~2+yl*s+y0) ; 

subs (x0=q0+7 , xl=ql+6*sqrt ( 15) /5 , y0=r0+27/125 , y l=r l+54*sqrt ( 15) /125 , 

y2=r2-43/5,s=t+a,a=-sqrt(15)/5,p) ; 
collect (%,t) ; 
simplify(°/J ; 
collect (%,t) 



show that the resulting polynomial is 

t^t 6 + A(q,r)(t) 
where the linear map A : C 5 — > P 5 is given by 

A(q,r)(t) 



(5) 



qit 5 + (g - v^gi) t 4 + \ 8 qi -^VlEq ^jt 3 + 

/28 12 r- \ 9 

( ygo-yVl5gi + r 2 Jt 2 + 

-32 /— 27 2 r— \ 

—— V15 g + -=-qi + n - -Vl5r 2 ) t + 
25 5 5 / 



39 
25 



go 



39 



— 1 i— 3 

. [9r . . 15 gi + r - -Vl5ri + -r 2 . 



It is easily verified that ^4(0, 0) = 0; hence the map A is indeed linear. Clearly, minimizing 
the abscissa of the polynomial (0) over [q, r] 6 C 5 is equivalent to the original problem 
(JTJ). Because the space of monic polynomials is not a linear space, it is convenient to 
introduce the notation 

70) = max{Re t : t n+l + w(t) = 0}, we P n , 

for the abscissa of t n+1 + w(t). We wish to establish that is a sharp local minimizer 
of the composition of the function 7 with the linear map A over [q, r] in the parameter 
space C 5 . 

To proceed further we need the notion of the adjoint map A* : P 5 — > C 5 , defined by 

'w{t), A{q,r)(t)) = (A*(w), [q, r] 



for all polynomials weP 5 and vectors [q, r] G C 5 , where the second inner product is the 
usual real inner product on C 5 and the first is a real inner product on P 5 , namely 

55 5 

3=0 3=0 ' 3=0 

It is easy to see that A* is given by 



A*(J2c3t j 

vj=0 
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2 5 Co - HVl5ci + fc 2 - § Vl5c 3 + c 4 



^Vl5co + f Cl -fVl5c 2 + 8c 3 
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15 c 4 + c 5 



-Iv^co + Ci 
\c - l\/l5ci + c 2 



(6) 



Following |4, and |3 j , we will establish that is a sharp local minimizer of the composition 
of 7 with the linear map A, which we denote 7 o A, by showing that 



G int 5(70 A) (0), 



(7) 



where d is the sub differential operator of variational analysis [3], [3 Chap. 8]. In order 
to do this we can use the nonsmooth chain rule [3J Lemma 4.4] 



a (70 a) = A*d 7 (o), 

as long as we verify the constraint qualification 

Af(A*) n <9°°7(0) = {0}, 



(8) 



(9) 



where M denotes null space and d°° denotes the horizon subdifferential operator [01 
Chap. 8]. This chain rule is valid because of the subdifferential regularity (HI Chap. 8] of 
the function 7 on P", established in 2j. The following formulas for the subdifferential and 
horizon subdifferential of 7 at were also established in but we follow the notation 
used in Theorem 3.3]: 



57(0) 
<9°°7(0) 



3=0 

n 

3=0 



n 



-, Re c n _i < 



0, Re c n _i < 



It follows from the latter formula that the constraint qualification (jOJ) holds if C5 = and 



implies c 



[ c 0; " " " j c 5\ 



0, a fact that is easily checked by observing 



that the 2 by 2 linear system 



fx/15 1 
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has only the trivial solution [03, C4] =0. Thus the nonsmooth chain rule 

5 



yields 



d( 1 oA) = {A*(J2^ 



: c 5 = -i Re c 4 < 0}. 
o 



The final step is to determine whether is in the interior of this subdifferential set. To 
check this, we need to solve the following linear system: set the right-hand side of © to 
as well as C5 



-|, which reduces to 



-fV15 



c 3 
c 4 





1 

6 



— -^VTE]. Since 



This linear system has a unique solution, namely [(%, C4] = [— ^ 
this satisfies the inequality Re C4 < 0, it follows that is in the subdifferential set 
and furthermore, since the inequality holds strictly, that every point near is in the 
subdifferential set, and therefore that (J7J) holds. This completes the proof of the theorem. 
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Step Response 
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Figure 2: Closed-loop step response. 

6 Time Response and Fragility 

In Figure El we graph the step response of the two mass-spring system fed back with 
controller (J3J, obtained with the following MATLAB [S] commands: 

P = tf (1, [1 2 0]); 

K = tf([-43/5 54*sqrt(15)/125 27/125] , [1 6*sqrt(15)/5 7]); 
T = feedback (P,K) ; 
step(T) ; 

We can see that the settling time is around 16 seconds. 

It is well known that multiple roots of polynomials are very sensitive to perturbations 
in the coefficients. In practice, this means that the closed-loop system will be fragile, 
or non- robust, with respect to uncertain data, implementation errors, or even rounding 
errors. For example, if instead of implementing the exact second-order controller ({5)1. we 
implement the nearby controller, obtained by keeping 5 significant digits, given by 

Xo = 7, xi = 4.6476, y = 0.2160, y x = 1.6731, y 2 = -8.6, 

then we obtain closed-loop poles at -0.9405, -0.8163 ± 0.1489 i, -0.7500 and -0.6622 ± 
0.0786 i, quite far from the single pole at —0.7746 assigned with the exact controller. 
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Figure 3: Possible zeros of the closed-loop polynomial p(x*,y*) under real perturbations 
to its coefficients with norm < 10~ 4 . 



This phemonenon can be studied graphically. In Figure EH we show 2 , in the gray region, all 
possible roots of polynomials that can be obtained by real perturbations to the polynomial 
p(x*, y*) given in (JJJ), where the norm of the vector of perturbations to the coefficients is 
no more than e = 1CT 4 . This region is sometimes called the real pseudozero set. 



7 Concluding Remarks 

In this note we have formulated the problem of maximizing the closed-loop asymptotic 
decay rate of a linear control system as a nonsmooth, nonconvex problem of polynomial 
abscissa minimization, focusing on the case of a benchmark two-mass-spring system. We 
derived a formula for a second-order controller with closed-loop poles clustered at a single 
point. Our main contribution is the use of recently developed techniques from nonsmooth 
variational analysis to prove local optimality of this controller. 

Motivated by this result, as well as the result in 0] on which it is based, very recent work 
PP using a completely different technique shows that the second-order controller described 
above is actually globally optimal. 

2 Thanks to S. Graillat, N. Higham and F. Tisseur for providing matlab scripts for the computation 
of real pseudozero sets of a polynomial. 
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Finally, it should be emphasized that asymptotic decay rate maximization is not, by itself, 
a practical objective. As shown graphically in Figure |HJ our locally optimal controller 
yields a closed-loop system which is sensitive to uncertainty and/or disturbance. In other 
words, lack of robustness is the price one has to pay to maximize the decay rate. In a 
typical control engineering system, a trade-off should be found between the asymptotic 
decay rate and other quantities, such as the complex or real stability radius, the complex or 
real pseudoabscissa (the maximum real part of the points in the complex or real pseudozero 
set), or H2 or performance measures. 
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